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The LEP2 experiments pose a serious naturalness problem for supersymmetric mod¬ 
els. The problem is stronger in gauge mediation than in supergravity models. Particu¬ 
lar scenarios, like electroweak baryogenesis or gauge mediation with light messengers, 
are strongly disfavoured. Searching a theoretical reason that naturally explains why 
supersymmetry has not been found poses strong requests on model building. If in¬ 
stead an unlikely (p ~ 5%) numerical accident has hidden supersymmetry to LEP2, 
we compute the naturalness distribution of values of allowed sparticle masses and 
supersymmetric loop effects. We find that b — » sy remains a very promising signal 
of minimal supersymmetry even if there is now a 20% (4%) probability that coloured 
particles are heavier than 1 TeV (3 TeV). We study how much other effects are 
expected to be detectable. 


1 Introduction 

Extensive searches of supersymmetric signals, done mostly 
but not only at LEP, have found no positive result so 
far [1, 2]. Nevertheless these results turn out to have in¬ 
teresting consequences, because the typical spectrum of ex¬ 
isting ‘conventional’ supersymmetric models contains some 
sparticle (a chargino, a neutralino or a slepton, depending 
on the model) with mass of few 10 GeV. How much should 
one worry about the fact that experimental bounds require 
these particles heavier than (80 H- 90) GeV? 

A first attempt to answer this question has been made 
in [3], using the fine tuning (‘FT’) parameter [4] as a quan¬ 
titative measure of naturalness. After including some im¬ 
portant one loop effects [5, 6] that alleviate the problem, 
in minimal supergravity the most recent bounds require a 
FT (defined as in [6]) greater than about 6 [6, 7], and a FT 
greater than about 20 to reach values of tan (3 < 2. 

Still it is not clear if having a FT larger than 6 is wor- 
ryingly unnatural. Setting an upper bound on the FT is a 
subjective choice. The FT can be large in some cases where 
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nothing of unnatural happens (i.e. dimensional transmuta¬ 
tion) [8]. Does a Z-boson mass with a strong dependence 
on the supersymmetric parameters indicate a problem for 
supersymmetry (SUSY)? The answer is yes: in this case 
the FT is related to naturalness, although in an indirect 
way [9]. Rather than discussing this kind of details, in 
this paper we approach the naturalness problem in a more 
direct way. 

When performing a generic random scanning of the 
SUSY parameters, the density distribution of the final re¬ 
sults has no particular meaning. For this reason the sam¬ 
pling spectra that turn out to be experimentally excluded 
are usually dropped. In order to clearly exhibit the natu¬ 
ralness problem, we make plots in which the sampling den¬ 
sity is proportional to the naturalness probability (we will 
discuss its relation with a correctly defined ‘fine tuning’ co¬ 
efficient), so that it is more probable to live in regions with 
higher density of sampling points. Plotting together the 
experimentally excluded and the few still allowed points 
we show in fig.s 2, 3 in a direct way how strong are the 
bounds on supersymmetry. We find that 

• The minimal FT can be as low as 6, but such a rela¬ 
tively low FT is atypical: the bounds on the chargino 
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and higgs masses exclude 95% of the MSSM parame¬ 
ter space with supersymmetry breaking mediated by 
‘minimal’ or ‘unified’ supergravity. 

• LHC experiments will explore 90% of the small re¬ 
maining part of the parameter space (we are assum¬ 
ing that it will be possible to discover coloured su¬ 
perparticles lighter than 2TeV). 

• The supersymmetric naturalness problem is more se¬ 
rious in gauge mediation models. 

Some regions of the parameter space are more problematic: 

• values of tan /3 lower than 2 (a range suggested by 
an infra-red fixed-point analysis and by the b/r uni¬ 
fication) ; 

• values of tan j3 bigger than 20 can naturally appear 
only in some particular situation, but often imply a 
too large effect in b —> sy. 

Some particular and interesting scenario appears strongly 
disfavoured 

• Gauge mediation models with low messenger mass 
M G m <, 10 7 GeV; 

or even too unnatural to be of physical interest 

• Baryogenesis at the electroweak scale allowed by a 
sufficiently light stop. 

These results are valid in ‘conventional’ supersymmetric 
models. It is maybe possible to avoid these conclusions 
by inventing appropriate models. However, if one believes 
in weak scale supersymmetry and thinks that the natu¬ 
ralness problem should receive a theoretical justification 
different from an unlikely numerical accident, one obtains 
strong constraints on model building. If instead super- 
symmetry has escaped experimental detection because a 
numerical accident makes the sparticles heavier than the 
natural expectation, we compute the naturalness distribu¬ 
tion probability for supersymmetric signals in each given 
model. There is a non negligible probability that an acci¬ 
dental cancellation stronger than the ‘minimal’ one makes 
the sparticles heavier than 1 TeV. 

In section 2 we outline and motivate the procedure that 
we adopt. In section 3 we show our results for the masses 
of supersymmetric particles. In section 4 we discuss the 
natural range of various interesting supersymmetric loop 
effects. Finally in section 5 we give our conclusions and we 
discuss some implication of the supersymmetric natural¬ 
ness problem for model building. In appendix A we collect 
the present direct experimental bounds on supersymmetry 
in a compact form. Since our formulas for supersymmet¬ 
ric effects in B-mixing correct previous ones and include 
recently computed QCD corrections, we list them in ap¬ 
pendix B. 


2 The procedure 

In this section we describe the sampling procedure we have 
used to make plots with density of points proportional to 
their naturalness probability. 

Within a given supersymmetric model (for example 
minimal supergravity) we extract random values for the 
dimensionless ratios of its various supersymmetry-breaking 
parameters. We leave free the overall supersymmetric mass 
scale, that we call msuSY ■ For each choice of the random 
parameters we compute v and tan /? by minimizing the po¬ 
tential: msuSY is thus fixed by imposing that v (or Mz) 


gets its experimental value. We can now compute all the 
masses and loop effects that we want to study; and com¬ 
pare them with the experimental bounds. 

This ‘Monte Carlo’ procedure computes how frequently 
numerical accidents can make the Z boson sufficiently lighter 
than the unobserved supersymmetric particles. As dis¬ 
cussed in the next subsection, in this way we obtain a sam¬ 
ple of supersymmetric spectra with density proportional to 
their naturalness probability. In subsection 2.2 we illus¬ 
trate these considerations with a simple example. 


2.1 Naturalness and fine tuning 


Our approach is related to previous work in this way: for 
any given value of the dimensionless ratios the natural¬ 
ness probability given by the procedure outlined before is 
inversely proportional to the ‘fine tuning-like’ parameter 
A defined in [9], that in the limit A 1 reduces to the 
original definition of the fine tuning parameter [4]. 

In a supersymmetric model the Z boson mass, Ml{p), 
can be computed as function of the parameters p of the 
model, by minimizing the potential. A Z boson much 
lighter than the supersymmetric particles is obtained for 
certain values p of the supersymmetric parameters, but is 
characteristic of only a small region Ap of the parameter 
space. Unless there is some reason that says that values in 
this particular range are more likely than the other possi¬ 
ble values, there is a small probability — proportional the 
size Ap of the small allowed range — that supersymme¬ 
try has escaped experimental detection due to some un¬ 
fortunate accident. This naturalness probability is thus 
inversely proportional to the ‘fine-tuning-like’ parameter 


A = 


P 

Ap 


(1) 


defined in [9]. As shown in [9], if Ap is small, we can 
approximate Ml(p) with a first order Taylor expansion 
around its experimental value, finding Ap « p/d[p] where 
d coincides with the original definition of fine-tuning [4] 


P dM\ 
M 2 Z dp 


( 2 ) 


This original definition of naturalness in terms of the loga¬ 
rithmic sensitivities of Mz has been criticized [8] as being 
too restrictive or unadequate at all. The definition (1) 
avoids all the criticisms [8]. 

To be more concrete, in the approximation that the Z 
boson mass is given by a sum of different supersymmetric 
contributions, the previous discussion reduce to saying that 
there is an unnaturally small probability p « A -1 that 
an accidental cancellations allows a single ‘contribution’ to 
Ml be A time bigger than their sum, Ml. 

Beyond studying how naturally Mz is produced by the 
higgs potential, we will study the natural expectation for 
various signals of supersymmetry. Our procedure automat¬ 
ically weights as more unnatural particular situations typi¬ 
cal only of restricted ranges of the parameter space that the 
ordinary ‘fine tuning’ parameter does not recognise as more 
unnatural (three examples of interesting more unnatural 
situations: a cancellation between too large chargino and 
higgs contributions to b —> sy\ an accidentally light stop; a 
resonance that makes the neutralino annihilation cross sec¬ 
tion atypically large). In the usual approach one needs to 
introduce extra fine-tuning coefficients (that for example 
measure how strong is the dependence of BR(6 —> sy), of 
the stop mass, of the neutralino relic density on the model 
parameters) to have a complete view of the situation. 
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Figure 1: The naturalness problem in a typical super- 
symmetric model. We plot the chargino mass in GeV 
as function of y/M 2 , that is the only free parameter 
of the model under consideration. Values of p marked 
in gray are unphysical, while light gray regions have 
too light sparticles. Only the small white vertical band 
remain experimentally acceptable. 



2.2 An example 

We now try to illustrate the previous discussions with a 
simple and characteristic example. We consider the ‘most 
minimal’ gauge mediation scenario with very heavy mes¬ 
sengers (one 5 © 5 multiplet of the unified gauge group 
SU(5) with mass Mgm = 10 15 GeV). ‘Most minimal’ means 
that we assume that the unknown mechanism that gener¬ 
ates the y-term does not give additional contributions to 
the other parameters of the higgs potential. Since the B 
terms vanishes at the messenger scale, we obtain moder¬ 
ately large values of tan /3 ~ 20. This model is a good 
example because its spectrum is not very different from a 
typical supergravity spectrum, and it is simple because it 
has only two free parameters: the overall scale of gauge me¬ 
diated soft terms and the y term. The condition of correct 
electroweak breaking fixes the overall mass scale, and only 
one parameter remains free. We choose it to be p = y/M 2 
(renormalized at low energy), where M 2 is the mass term 
of the SU(2)i gaugino. In this model we can assume that y 
and M 2 are real and positive without loss of generality. In 
figure 2.2 we plot, as a function of p, the lightest chargino 
mass in GeV*. The gray regions are excluded because the 
electroweak gauge symmetry cannot be properly broken if 
y/M 2 is too small or too large. Values of p shaded in light 
gray are excluded because some supersymmetric particle is 
too light (in this example the chargino gives the strongest 
bound). The chargino is heavier than its LEP2 bound only 
in a small range of the parameter space at p ss 2.3 (left 
unshaded in fig. 1) close to the region where EW symmetry 
breaking is not possible because y is too large. The fact 

I We must mention one uninteresting technical detail since 
figure 1 could be misleading about it. We are studying how 
frequently numerical accidents can make the Z boson lighter 
than the sparticles. We are not studying how frequently nu¬ 
merical accidents can make the sparticles heavier than the Z 
bosons. When including loop effects the two questions are not 
equivalent. The first option, discussed in this paper because 
of physical interest, gives sparticles naturally heavier than the 
second option. 


that the allowed regions are very small and atypical is a 
naturalness problem for supersymmetry. Roughly 95% of 
the possible values of p are excluded by the chargino mass 
bound. 

In this example we have fixed M t pole = 175 GeV, and we 
have taken into account the various one-loop corrections to 
the effective potential (that double the size of the allowed 
parameter space). In our example we have used a small 
value of the top quark Yukawa coupling at the unification 
scale, At (VIgut) ~ 0.45, not close to its infrared fixed point 
but compatible with the measured top mass. In this way 
the large RGE corrections to the higgs masses are mini¬ 
mized. This directly alleviates the naturalness problem; 
moreover, since the B term is only generated radiatively, 
we also get a moderately large value of tan/3 ~ 20 that 
indirectly further alleviates the naturalness problem. As a 
consequence the fine tuning corresponding to this example 
is A = 6, the lowest possible value in unified supergravity 
and gauge mediation models. In the next section we will 
study the naturalness problem in these motivated models. 
To be conservative, up to section 4 we will use only accel¬ 
erator bounds, and we will not impose the constraint on 
the b —* sy decay (and on other loop effects). 

3 The supersymmetric natural¬ 
ness problem 

In this section we discuss how serious is the supersymmet¬ 
ric naturalness problem in the various different motivated 
scenarios of supersymmetry breaking. We will not consider 
models with extra fields at low energy beyond the minimal 
ones present in the MSSM and we will assume that ‘matter 
parity’ (equivalent to R parity) is conserved. 

In all this section we will consider as excluded only 
those spectra that violate the experimental bounds coming 
from direct searches at accelerators listed in appendix A. 
We do not impose cosmological bounds (the addition of 
tiny R parity violating couplings allows to remove even¬ 
tual problems of dark matter overabundance or of nucle¬ 
osynthesis destruction); we do not impose that the physical 
vacuum be the only (or the deepest) one [10, 11] (depend¬ 
ing on cosmology the presence of extra unphysical minima 
can be or cannot be dangerous — within conventional cos¬ 
mology (i.e. sufficiently hot universe, inflation) the most 
frequent unphysical minima seem not dangerous [11]); we 
do not impose b/r unification or closeness of At to its infra¬ 
red fixed point (these appealing assumptions give prob¬ 
lems [7, 12] that can be alleviated by modifying the the¬ 
ory at very high energy or by loop corrections at the elec¬ 
troweak scale); and in this section we also do not impose 
any indirect bound (because we want to be conservative 
and include only completely safe bounds — the indirect 
constraint from BR(£? —> X s y) would exclude 20% of the 
otherwise allowed points). We are thus excluding spectra 
that are really excluded. 


3.1 Minimal supergravity 

“Minimal supergravity” assumes that all the sfermion mas¬ 
ses have a common value mo, all the three gaugino masses 
have a common value Ms, and that all the A-terms have 
a common value Aq at the unification scale Mgut ~ 2 • 
10 16 GeV. The parameters yo and Bo are free. They con¬ 
tribute to the mass terms of the higgs doublets h u and h d 
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Figure 2: Scatter plot with sampling density proportional to the naturalness probability. The area shaded in 
light gray (dark gray) in fig.s 2a,b correspond to regions of each plane excluded at LEP2 (at LEP1), while the 
area shaded in light gray in fig. 2c has been excluded at Tevatron. The dark gray (black) points correspond to 
sampling spectra excluded at LEP2 (at LEP1). Only the light gray points in the unshaded area satisfy all the 
accelerator bounds. Points with unbroken electroweak symmetry are not included in this analysis. 


in the following way: 

(mfu + |/io|)|h u | 2 + (m 2 a + |^ol)|h d | 2 + {goBoh u h d + h.c.). 

As explained in the previous section we randomly fix the 
dimensionless ratios of the parameters mo, Ms, Aq, Bo 
and go and we fix the overall supersymmetric mass scale 
“msuSY” from the minimization condition of the MSSM 
potential. More precisely we scan the parameters within 
the following ranges 


mo 

= (g 4- 3)msusY 

(3a) 

\go\,M 5 

= (- 4- 3)msusY 

(3b) 

Ao, Bo 

= (—3 4- 3)m 0 

(3c) 


The samplings in (3a), (3b) are done with flat density 
in logarithmic scale. We think to have chosen a reasonable 
restriction on the parameter space. We could make the nat¬ 
uralness problem apparently more dramatic by restricting 
the dimensionless ratios to a narrow region that does not 
include some significant part of the experimentally allowed 
region, or by extending the range to include larger values 
that produce a larger spread in the spectrum so that it is 
more difficult to satisfy all the experimental bounds. The 
only way to make the naturalness problem less dramatic is 
by imposing appropriate correlations among the parame¬ 
ters, but this makes sense only if a theoretical justification 
can be found for these relations.We will comment about 
this possibility in the conclusions. An alternative scanning 
procedure that does not restrict at all the parameter space 

mo, |mo|, Ms, Bo, |Ao| = (0 4- l)msusY- 

gives the same final results as in (3). 

We now exhibit the results of this analysis in a se¬ 
ries of figures. In fig.s 2 we show some scatter plots with 
sampling density proportional to the naturalness proba¬ 
bility. The sampling points that give spectra excluded at 
LEP1 are drawn in black, the points excluded at LEP2 in 
medium gray, and the still allowed points are drawn in light 
gray. The present bounds are listed in appendix A. The 


pre-LEP2 bounds approximately consist in requiring that 
all charged and coloured particles be heavier than Mz / 2 . 
Fig. 2a shows the correlation between the masses of the 
lightest chargino, m x , and of the lightest higgs, This 
plot shows that the experimental bounds on m x and mh 
are the only important ones in minimal supergravity (if 
mo M 5 also the bound on the mass of right-handed 
sleptons becomes relevant). The bound on the chargino 
mass is more important than the one on the higgs mass: 
even omitting the bound on the higgs mass the number 
of allowed points would not be significantly increased (the 
MSSM predicts a light higgs; but this prediction can be 
relaxed by adding of a singlet field to the MSSM spec¬ 
trum). The experimental bounds on supersymmetry are 
thus very significant (to appreciate this fact one must no¬ 
tice that the allowed points have small density). This fact 
is maybe illustrated in a more explicit way in fig. 2 b, where 
we show the correlation between the masses of the lightest 
chargino, M x , and of the lightest neutralino, Mn■ We see 
that in the few points where LEP2 bound on the chargino 
mass is satished, the two masses Mn ~ Mi and M x ss M 2 
are strongly correlated because the g parameter is so large 
that the SU(2)i-breaking terms in the gaugino mass ma¬ 
trices become irrelevant. We clearly appreciate how strong 
has been the improvement done at LEP2. In fig. 2c we 
show an analogous plot in the {Mo,mf) plane, often used 
to show bounds from hadronic accelerators. The bounds 
from LEP2 experiments together with the assumption of 
mass universality at the unification scale gives constraints 
on the mass of coloured particles stronger than the direct 
bounds from Tevatron experiments 

In fig.s 3 we show the same kind of results using a 
different format. We show the ‘naturalness distribution 
probability’ for the masses of various representative super- 
symmetric particles and for tan j3. The allowed points have 
been drawn in dark gray; those ones excluded at LEP2 
(LEP1) in in medium gray (in light gray). As a conse¬ 
quence of the naturalness problem the allowed spectrum 
is confined in fig.s 3to the small tails in the upper part 
of the unconstrained probability distributions. The fact 
that LEP1 experiments have not found a chargino lighter 
than Mz /2 excluded about 70% of the unified supergrav- 
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Figure 3: Naturalness distribution of sparticle masses in minimal supergravity. Allowed spectra contribute only 
with the small tails in dark gray. The remaining 95% of the various bell-shaped distributions is given by points 
excluded at LEP2 (in medium gray) or at LEP1 (in light gray). On the vertical axes on each plot the particle 
masses in GeV are reported (tan /? in the first plot). With ‘squark’ and ‘slepton’ me mean the lightest squark 
and slepton excluding the third generation ones, that have weaker accelerator bounds. 


ity parameter space. In the MSSM with soft terms me¬ 
diated by ‘minimal supergravity’ the present experimental 
bounds exclude 95% of the parameter space with broken 
electroweak symmetry (97% if we had neglected one-loop 
corrections to the potential). 

3.2 Non minimal supergravity 

Some of the assumptions of the ‘minimal supergravity sce¬ 
nario’ allow to reduce the number of parameters, but do 
not have a solid theoretical justification. ‘Minimal super¬ 
gravity’ at the Planck (or string) scale can be justified: 
but even in this case the mass spectrum at the unification 
scale is expected to be significantly different from the min¬ 
imal one due to renormalization effects [13]. If the theory 
above the gauge unification scale is a unified theory, these 
renormalization effects induce new flavour and CP violat¬ 
ing processes [14, 15], that we will discuss in section 4.2. In 
this section we study the naturalness problem in ‘unified 
supergravity’, that we consider an interesting and predic¬ 
tive scenario. More precisely, beyond performing the scan¬ 
ning (3), we also allow the higgs mass parameters at the 
unification scale to vary in the range 

m h u,m h a = (i 4- 3)m 0 

where mo is now the mass of the sfermions contained in the 
IO 3 = ((In, bif), tn, tr) multiplet of SU(5). We can assume 
that all remaining sfermions have the same mass mo, since 
they do not play any significant role in the determination 
of M§ (unless they are very heavy; we do not consider this 
case in this section). 

From the point of view of the naturalness problem there 
is no significant difference between minimal and unified 
supergravity: an unnatural cancellation remains necessary 
even if there are more parameters. 


From the point of view of phenomenology, maybe the 
most interesting new possibility is that a mainly right handed 
stop can ‘accidentally’ become significantly lighter than the 
other squarks. Unless the top ^4-term is very large, this ac¬ 
cidental cancellation is possible only if [17] 

mfo 3 /wh u ~ (0.6 4- 0.8) and M 2 <) 0.2|m/j u | 

(all parameters renormalized at Mgut) — a region of the 
parameter space where the small mass M 2 term for the 
chargino makes the naturalness problem stronger. The 
large Yukawa coupling of a light stop can generate various 
interesting loop effects ( b —> s'y, A ms, ek) and make the 
electroweak phase transition sufficiently first-order so that 
a complex g term can induce baryogenesis [16]. However 
this possibility is severely limited by naturalness considera¬ 
tions [17]: it requires not only that one numerical accident 
makes the Z boson sufficiently lighter than the unobserved 
chargino, but also that a second independent numerical 
accident gives a stop lighter than the other squarks. The 
fact that in our analysis this combination occurs very rarely 
(p ~ 10 -3 ) means that this nice possibility is very unnat¬ 
ural. Furthermore a stop lighter than 200 GeV is strongly 
correlated with a too large supersymmetric correction to 
the BR(B —> Xs'y) decay. 

In fig. 4 we show the naturalness distribution of sparti¬ 
cle masses omitting all the experimentally excluded spec¬ 
tra, i.e. under the assumption that supersymmetry has es¬ 
caped detection because a numerical accident makes the 
sparticle too heavy for LEP2. In the more realistic ‘unified 
supergravity’ scenario, these distributions correspond to 
the ‘allowed tails’ of the full distributions plotted in fig.s 3. 
Since the allowed spectra come from a small region of the 
parameter space, their naturalness distribution probabil¬ 
ity has a less significant dependence on the choice of the 
scanning procedure. 
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It is possible to define the most likely range of values of 
the masses of the various sparticles, for example by exclud¬ 
ing the first 10% and the last 10% of the various distribu¬ 
tions. These (10 4- 90)% mass ranges are shown in table 1. 
We see that, once one accepts the presence of a numerical 
accident, it is not extremely unlikely that it is so strong 
that coloured particles are heavier than few TeV (even for 
so heavy sparticles, a stable neutralino is not necessarily 
dangerous for cosmology). Nevertheless there is still a very 
good probability that supersymmetry can be detected at 
LHC . In section 4 we will use these naturalness probability 
distributions for sparticle masses to compute the natural 
values of various interesting loop effects mediated by the 
sparticles, that could be discovered before LHC. 

Before concluding this section, we recall that the bounds 
on the gluino mass from Tevatron experiments are not com¬ 
petitive with the LEP2 bounds on the chargino mass, if 
gaugino mass universality is assumed. Consequently it is 
possible to alleviate the naturalness problem of supersym¬ 
metry by abandoning gaugino mass universality and mak¬ 
ing the gluino pole mass as light as possible, M 3 ~ (200 4- 
250) GeV (in supergravity this requires a small gluino mass, 
Ms(Mgut) ~ 100 GeV, at the unification scale; an analo¬ 
gous possibility in gauge-mediation requires a non unified 
spectrum of messengers). In this case the sparticle spec¬ 
trum is very different from the one typical of all conven¬ 
tional models; consequently the loop effects mediated by 
coloured sparticles are larger. Even if in this case the min¬ 
imal FT can be reduced even down to 0(1 ) values, we do 
not find this solution completely satisfactory. If we treat 
all the gaugino and sfermion masses as free parameters of 
order Mz, the SM Z boson mass is still one (combination) 
of them and still there is no reason, different from an un¬ 
welcome accident, that explains why Mz is roughly the 
smallest of all the ~ 10 charged and coloured soft masses. 
Anyway, Tevatron can concretely explore this possibility in 
the next years. 

3.3 Gauge mediation 

‘Gauge mediation’ models contain some charged ‘messen¬ 
ger’ superfields with some unknown mass Mgm = (10 5 4- 
10 15 )GeV directly coupled to a gauge singlet field with 
supersymmetry-breaking vacuum expectation value. The 
‘messengers’ mediate supersymmetry-breaking terms to the 
MSSM sparticles that feel gauge interactions [18, 19]. The 


supersymmetric 

particle 

mass range 
(in GeV) 

lighest neutralino 

55 4- 250 

lighest chargino 

110 4- 500 

gluino 

400 4- 1700 

slepton 

105 4- 600 

squark 

400 4- 1700 

stop 

250 4- 1200 

charged higgs 

300 4- 1200 


Table 1: (10 4-90)% naturalness ranges for the masses 
of various supersymmetric particle, in unified su¬ 
pergravity, assuming that the naturalness problem is 
caused by an accidental cancellation. 


spectrum of the supersymmetric particles is thus mainly 
determined by their gauge charges. More precisely, in a 
large class of minimal models (from the first toy ones to 
the more elaborated recent ones) the prediction for the soft 
terms, renormalized at the messenger mass Mgm, can be 
conveniently parametrized as 

Mi (M gm ) = Q ^ Gm) M 0 , (4a) 

47T 

m 2 R (M GM ) ~ g ■ CjjM^Mgm), (4b) 

where m R are the soft mass terms for the fields 
R = Q,u R ,d R ,e R ,L,h u ,h , 

and the various quadratic Casimir coefficients d R are fisted, 
for example, in ref. [9]. Here Mo is an overall mass scale 
and g parametrizes the different minimal models. For ex¬ 
ample g = (ns + 3nio) -1/,:2 < 1 in models where a single 
gauge singlet couples supersymmetry breaking to ns copies 
of messenger fields in the 5© 5 representation of SU(5) and 
to n 10 copies in the 10 © 10 representation [18, 19]. Val¬ 
ues of g bigger than one are possible if more than one 
supersymmetry-breaking singlet is present [18, 19]. If the 
messengers are as light as possible the sparticle spectrum 
become a bit different than the one in (4) (but we will 
argue in the following that this case is very unnatural). 
Gauge-mediation models have the problem that gauge in¬ 
teractions alone cannot mediate the ‘/t-term’, as well as 
the corresponding ‘B ■ /x-term’, since these terms break a 
Peccei-Quinn symmetry. The unknown physics required 
to solve this problem may easily give rise to unknown non 
minimal contributions to the soft terms in the Higgs sec¬ 
tor [19], but this lack of predictivity does not prevent the 
study of naturalness. 

Like in supergravity models, the one loop corrections to 
the potential are very important. In gauge mediation mod¬ 
els the minimal fine tuning is higher than in supergravity 
models [9, 20] because gauge mediation generates a right 
handed selectron mass significantly smaller than the higgs 
mass term which sets the scale of electroweak symmetry 
breaking. This effect is more pronounced for intermediate 
values of the mediation scale, Mgm ~ 10 8 GeV. Higher val¬ 
ues of the mediation scale give spectra of sparticles more 
similar to supergravity case, that have a less strong nat¬ 
uralness problem (but if Mgm <; 10 12 GeV it is necessary 
to complicate the theory to avoid destruction of nucleosyn¬ 
thesis [19]). If the mediation scale is as light as possible the 
RGE effects between Mgm and Mz are smaller, reducing 
the fine-tuning. However for low values of the messenger 
mass, Mgm ^ 10'~ 8 GeV, the neutralino decays within the 
detector, so that LEP2 experiments give now a very strong 
bound on its mass: mjv > 91 GeV (if 77 ^ 0.5 the experi¬ 
mental bounds can be less stringent because a slepton could 
be lighter than the neutralinos). The consequent strong 
naturalness problem makes the light messenger scenario 
not attractive; Values of the gauge mediation scale higher 
than Mgm 10 9 GeV are instead less attractive from the 
point of view of cosmology [21]. 

In fig. 5 we show contour plots of the fine tuning re¬ 
quired by gauge mediation models, in the plane (Mgm,??) 
for tan (3 = 2.5 and M t pole = 175 GeV (this corresponds 
to A((Mgut) ~ 0.5, taking into account threshold correc¬ 
tions to At; the FT is somewhat higher if one uses larger 
allowed values of At(MGUT)), where Mgm is the messen¬ 
ger mass, while g parametrizes the different models, as de¬ 
fined in (4). The fine tuning is computed according to the 
definition in [9], including one loop effects and all recent 
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Figure 4: Naturalness distribution of some illustrative supersymmetric particle masses (on the horizontal axis) 
in unified supergravity, under the hypothesis that supersymmetry has not be found at LEP2 due to a numerical 
accident. The three continuous lines are the three gauginos (Mn, M x and Mg from left to right). The thin 
dashed line is the lightest slepton and the thick dashed lines are the lightest stop (left) and squark (right). The 
dotted lines are the light and charged higgs. The distribution of tan /3 is not shown because similar to the one 
in fig. 3. 


bounds. In these figures we have shaded in medium gray 
the regions at high values of the messenger mass where the 
bound on the chargino mass is the strongest one and the 
fine-tuning is not higher than in supergravity models. We 
have shaded in dark gray the regions at small values of the 
messenger mass where the very strong new bound on the 
neutralino mass makes the model unattractive. We have 
not coloured the remaining region where the bound on the 
right handed slepton masses is the strongest one and makes 
the fine tuning higher than in supergravity models. Small 
(?7<[0.4) and big (g > 1) values of g (below and above 
the dashed lines) are only allowed in models where more 
than one singlet field couples supersymmetry breaking to 
the messenger fields. The fine tuning strongly increases 
for smaller values of tan /3, and becomes a bit lower for 
higher values of tan j3. The choice of parameters used in 
the example shown in fig. 1 corresponds to one of the most 
natural cases. 

As indicated by fig. 5, extremely light messengers (Mgm ~ 
10 TeV) could give rise to a more natural sparticle spec¬ 
trum (with detectably non-unified gaugino masses [22]); 
we have not studied this possibility because NLO correc¬ 
tions [22, 23], that depend on unknown couplings between 
messengers, become relevant in this limiting case. 

4 Supersymmetric loop effects 

Assuming that the supersymmetric naturalness problem 
is caused by a numerical accident, we now study the natu¬ 
ral values of various supersymmetric loop effects, hopefully 
detectable in experiments at energies below the supersym¬ 
metric scale. As discussed before, since the allowed param¬ 
eter space is very small, we can safely compute naturalness 
distributions in any given model. As before we concentrate 
our analysis on unified supergravity and, for simplicity, we 
assume that the sfermions of each given generation have 
a common soft mass and a common A term at the the 
unification scale. In a first subsection we assume a com¬ 
plete degeneration of all the sfermions at the the unification 
scale. In this case the CKM matrix is the only source of 
flavour and CP violation (possible complex phases in the 
supersymmetric parameters — not discussed here — would 
manifest at low energy mainly as electric dipoles). The su¬ 


persymmetric corrections to loop effects already present in 
the SM are studied in this first subsection. 

In the second subsection we will assume (as suggested 
by unification of gauge couplings and by the heaviness of 
the top quark [14, 15]) that the sfermions of third gen¬ 
eration are non degenerate with the other ones. In this 
case the MSSM Lagrangian contains new terms that vi¬ 
olate leptonic flavour, hadronic flavour and CP and that 
can manifest in a variety of ways. These effects depend in 
a significant way not only on the masses of the sparticles, 
but also on unknown parameters not constrained by nat¬ 
uralness. We will thus compute them only in particular 
motivated models. 

4.1 Minimal effects 

In cases where electroweak loop corrections give observ¬ 
able effects to some measurable quantities supersymme¬ 
try can give additional loop corrections comparable to the 
SM ones. None of these corrections have been seen in the 
electroweak precision measurements done mostly at LEP. 
However more powerful tests of this kind will be provided 
by more precise measurement and more precise SM pre¬ 
dictions of some ‘rare’ processes in B physics (b —> sj, 
b —► s£ + £~, A ms), K physics ( ek , K —► Ttvv) and of 
the anomalous magnetic moment of the muon, = (g — 
2)^/2*. Notwithstanding the stringent constraints on the 

We have carefully computed the supersymmetric corrections 
to b — ► S 7 (assuming an experimental cutoff > 70 %E™ ax on 
the photon energy) including all relevant NLO QCD corrections 
in this way: The SM value have been computed as in [24] . The 
charged Higgs contribution is computed including only the rel¬ 
evant NLO terms as in the first reference in [25] (the other two 
references in [25] include all the remaining terms requested by 
a formal NLO expansion, that at most affect BR (B —> Xg'y) at 
the 1% level). The NLO corrections to the chargino/stop con¬ 
tribution are taken from [26] , again omitting the negligible NLO 
corrections to the b —► sg chromo-magnetic penguin. Unfortu¬ 
nately the approximation used in [26] is not entirely satisfactory: 
the chargino/squark corrections to the b —► S 7 rate can be large 
even without a lighter stop with small mixing, as assumed in [26] 
to simplify the very complex computation. Finally, we have in¬ 
cluded the threshold part of the two-loop O(X^) corrections that 
does not decouple in the limit of heavy supersymmetric parti¬ 
cles. We have computed using the formulae given in [27] . The 
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Figure 5: Contour plot of the fine-tuning parameter A in gauge mediation models in the plane f or 

p < 0 (left) and p > 0 (right). In the (dark gray, white, light gray) area at the (left, center, right) of each 
picture the strongest experimental constraint is the one on the mass of the lightest (neutralino, right handed 
slepton, chargino). Values of r) in the gray area below the lower dashed lines and above the upper dashed lines 
are not allowed in models with only one SUSY-breaking singlet. 


sparticle masses, in some particular region of the unified su¬ 
pergravity parameter space, all these supersymmetric cor¬ 
rections can still be significant [29]. We study how natural 
are these regions where the effects are maximal. 

In fig. 6 we show the naturalness distribution proba¬ 
bilities for these effects. We have drawn in light gray the 
contributions from spectra where one (or more) supersym¬ 
metric loop effect is too large; in medium gray the contribu¬ 
tions from spectra where one (or more) effect is detectable 
in planned experiments, and in black the contributions 
from spectra for which all effects are too small. Requir¬ 
ing a ~ 95% confidence level for considering excluded or 
discovered an effect, we divide the possible supersymmetric 
effects into ‘excluded’, ‘discoverable’, or ‘too small’ accord¬ 
ing the following criteria. We consider allowed a value of 
R = BR(B —> X s 7)mssm/BR(B —► X s 7)sm (where the 
MSSM value includes SM and sparticle contributions) in 
the range 0.6 < R < 1.5, and we consider discoverable a 
correction to R larger than 20%. The present experimen¬ 
tal uncertainty on 10 10 a M is ±84 [1] and will be reduced 
maybe even below the ±4 level [30]. It seems possible to 
reduce the present QCD uncertainty in the SM prediction 
to the same level [31]. In our plot we consider detectable 
a supersymmetric correction to a M larger than 10 -9 . It is 
easy to imagine how fig. 6 would change with more or less 
prudent estimates. 

The results are the following. Fig. 6 shows that b —> 
S 7 is a very promising candidate for a discoverable ef¬ 
fect. More precisely supersymmetry gives a correction to 
BR (B —> X s y) larger than 20% in about 40% of the al¬ 
lowed sampling spectra. A supersymmetric correction to 
the b —> sy magnetic penguin also manifests as a distortion 
of the spectrum of the leptons in the decay b —> s£ + £~. A 
detectable supersymmetric effect in the g — 2 of the muon 

supersymmetric correction to A mg have been taken from [28]. 


is less likely but not impossible. About 10% of the sampled 
spectra are accompanied by a discoverably large effect in 
g— 2 without a too large correction to b —> sy. On the con¬ 
trary a detectable (i.e. ]> 30%) supersymmetric correction 
to B mixing (shown in fig. 6c) and to K mixing (not shown 
because linearly correlated to the effect in B mixing) can 
only be obtained for values of the parameters [29] strongly 
disfavoured by our naturalness considerations. The same 
conclusion holds for K —» nui/ decays. 

An enhancement of the effects is possible in two par¬ 
ticular situations: if tan /3 is large, or if a stop state is 
lighter than ~ 200 GeV. Both these situation can be real¬ 
ized — but only for particular values of the parameters — 
in the ‘unified supergravity’ scenario in which we are doing 
our computations. As discussed in the previous section, a 
so light stop is decidedly not a natural expectation. The 
possibility of a large tan /3 is instead a weak aspect of our 
analysis. If the scalar masses are larger than all the other 
soft terms ( A , B, p and gaugino masses) a large tan /3 is 
naturally obtained [32]. In our scanning of the parameter 
space we have preferred to assume that the A terms are 
of the same order of the scalar masses, and we have thus 
not covered this possibility. We do not explore the pos¬ 
sibility of large tan (3 in this article because a large tan /3 
would enhance the one loop effects that are already more 
(b —> sy) or less (a M ) promising, but not the effects that 
seem uninteresting. 

4.2 New supersymmetric effects 

The mass matrices of the sfermions can contain new sources 
of flavour and CP violation, both in the hadronic and in 
the leptonic sector, that manifest themselves in processes 
either absent (like g —> ey) or extremely small (like the 
electron and neutron electric dipoles) in the SM. 
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Figure 6: Naturalness distribution of three possibly interesting ‘minimal’ supersymmetric effects (BR(.B —> 
* s 7)mssm/BR(B —> ^s7)sm ^ fig- 6a, anomalous magnetic moment of the muon a M in fig. 6b and 

Am^ USY /AmB in fig. 6c) in unified supergravity with degenerate sfermions. We have plotted in light gray 

the contributions from spectra where at least one of these loop effects is too large, in medium gray the ones 

where at least one of these loop effects is detectable, and in black the ones where all effects seem too small. Our 

boundaries between excluded/detectable/too small effects are represented by the horizontal lines. 


This possibility is often discussed (see e.g. [33]) in the 
‘mass insertion’ approximation [14], where the sfermion 
mass matrices are proportional to the unit matrix, plus 
small (and unknown) off diagonal terms. There is however 
no phenomenological constraint that forces the sfermions of 
third generation to be degenerate with the corresponding 
ones of the first two generations. Indeed, even with a max¬ 
imal 12/3 splitting, fermion/sfermion mixing angles 
as large as the CKM ones do not necessarily produce too 
large effects. Thus we will allow the masses of sfermions of 
third generation, m 3 j, to be different from the other ones, 
and parametrize this 12/3 non degeneration introducing a 
parameter pf. 

= m^/pf at the unification scale. (5) 

Rather than present a general parametrization, we now 
prefer to restrict our analysis to the case that we consider 
more strongly motivated [15]: order one 12/3 splitting (i.e. 
Pf ~ 1/2) and 12/3 mixing angles of the order of the CKM 
ones (i.e. Vfj ~ Vckm)*- For simplicity we continue to 
assume that all the sfermions of each generation have a 
common mass the unification scale (i.e. = mi, so that 

Pf = p) and we assume that no new effect of this kind 
comes from the A terms. 

* This scenario is motivated by the following considerations. 
The largeness of the top quark Yukawa coupling, At, suggests 
that the unknown flavour physics distinguishes the third gen¬ 
eration from the other ones. This is for example the case of a 
U(2) flavour symmetry [35]. A stronger motivation for pf yf 1 
comes from unification [14, 15]: the running of the soft terms in 
a unified theory gives pf < 1 even ii pf = 1 at tree level, due 
to the large value of the unified top quark Yukawa coupling. 
If A / ( A/q i;T) 7 1 this minimal effect is always very large [34], 
However these values of ARMqut) close to its IR fixed point 
accommodate the measured top mass only for tan 0<>2 — a 
range now disfavoured by the higgs mass bound together with 
naturalness consideration. At larger tan /l > 2 a top mass in 
the measured range requires ARMqut) = 0.35 -4- 0.55: for this 
smaller value the RGE running of the unified soft terms gives 
p ~ 0.8. More precisely, depending on the details of the model 
the effect in p can be large (p = 0.5) or very small ( p = 0.95). 


In the limit pf —> 0 only the third generation sparticles 
have mass around the electroweak scale. If also the third 
generation sleptons have few TeV masses, we encounter 
the scenario named as ‘effective supersymmetry’ in [36]. 
In the opposit limit, pf —> 1, we reduce to the ‘mass in¬ 
sertion’ approximation [14]. In both these cases the loop 
functions relevant for the various effects reduce to particu¬ 
lar limits. Since we consider non degenerate sfermions, we 
cannot use the ‘mass insertion’ parametrization (it is triv¬ 
ial to generalize it; but it becomes cumbersome since in 
some cases the dominant effects come from diagrams with 
two or three mass insertions). We choose pf = 0.5 (i.e. all 
the sfermions of third generation are significantly lighter 
than the other ones, a maybe too optimistic assumption) 
and we allow the fermion/sfermion mixing angles to vary 
in the range 

\Ve L r L \,\Ve RfR \,\V dL - b J,\V dRiR \ = (± - 3)[V*„| (6a) 

\V» L f L \,\V, R r R \,\V SL ~ bL \,\V aRbR \ = (i 4-3)[Vt s | (6b) 

(all angles are renormalized at the unification scale; we 
assume that at the weak scale Vts = 0.04 and Vts = 0.01) 
with complex phases of order one. Mixing angles in the up- 
quark sector are less motivated and thus less controllable. 
Their possible contribution to the neutron EDM will not 
be discussed here. 

Studying this case is a good starting point for under¬ 
standing what happens in similar cases. At the end of 
this subsection we will comment on how our results change 
if one of our simplifying but questionable assumptions is 
abandoned. 

Having made specific (but motivated) assumptions, we 
can now compute the resulting supersymmetric effects^ us¬ 
ing our naturalness distribution of sparticle masses. The 

'The supersymmetric effects are computed as follows. We 
take the expressions of the leptonic observables from [34] and 
the hadronic ones from [38] (where some of them are only given 
in symplifying limits); the supersymmetric CP asymmetry in 
b —> sss is taken from [39]. We have assumed, as in [39], that 
the largely unknown BR( B,i —> rf>Kg) is 10 -5 — a value con- 
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Figure 7: In fig.s 7a,b we show the ‘best fit’ values of the plane (p,r\). In fig. 7a we include all data in the fit, 
while in fig. 7b we omit Ek- In fig. 7c we include all data except Ex, and we study the compatibility between 
theory and experiments of each given value of p and g. In all cases the contour levels correspond to 68%, 95% 
and 99% C.L. 


results are shown in fig.s 8. The labels below each plot 
indicate the content of the plot; a longer description is 
written in the caption. We define here precisely the con¬ 
tent of the two plots at the right. The new supersymmetric 
effects that we are studying do not affect in a significant 
way BR(fe —> sy) but can modify in a detectable way its 
chiral structure. In the up-right plot we have plotted the 
supersymmetric contribution to the direct CP asymmetry 
l4u S7 | [41] defined by 


F(Bj — > X s y) — F(B d —> Xsy) 

dir r(B d x s y) + r(Ba - Xrt) B7>0 . 7Br * ■ 

In the down-right plot we have plotted the supersymmet¬ 
ric contribution to the mixing induced CP asymmetry [42] 
KTx S7 | defined by 


F(B 3 (t) -> My) - r (B d (t) My) 
T(B d (t) - My) + T (B d (t) - My) 


±-4mTx S7 sm([Am B Jt) 


where ± is the CP eigenvalue of the CP eigenstate M. If 
the Wolfenstein parameters p and g are given by the fit in 
fig. 7a, in the SM « +0.5% and M^x* 7 « 5%. It 

is however difficult to find particles M that allow a precise 
measurement of 4l^ S7 [42], 

We see that a non zero BR(/i —> ey) (proportional to 
BR(/i —> eee) and to the similar effect of p, —> e con¬ 
version [34], and strongly correlated with d e ) is the most 
promising candidate for a detectable effect. For a precise 
interpretation of the results we first discuss how large the 
various effects are now allowed to be, and the future ex¬ 
perimental prospects. 

sistent with the most reliable phenomenological estimate [40]. 
To compute the CP asymmetries in b — > sy we use the general 
formulae in [41, 42], In the computation of A mg and sk we add 
some important corrections with respect to previous analyses: 
we include the recently correctly computed QCD corrections [43] 
and the very recent lattice values of the matrix elements of the 
A B, AS = 2 supersymmetric effective operators [44, 45]. More¬ 
over we correct an error in the supersymmetric Wilson coefficient 
for A mg in [38, 39] (for a concidence the error only causes an 
irrelevant flip of the sign of the effect in the semi-realistic limit 
M 3 = mf). See appendix B for the details. 


Allowed supersymmetric effects 

Many of the experimental bounds are the same as the ones 
quoted in [34, 38]. Only the bound on the p —> ey decay 
and to p —> e conversion have been slightly improved by the 
Mega and SindrumII experiments. The present bounds 
are BR(u —> ey) < 3.8 1CF 11 [47] and CR(/rTi —> eTi) < 
6.1 lCT i3 [48], 

It is less clear how large the various present experi¬ 
mental data allow to be a supersymmetric correction to ek 
(i.e. to CP violation in KK mixing). The detailed study of 
this question is interesting because low energy QCD effects 
enhance the supersymmetric contribution to Ek, that for 
this reason becomes one of the more promising hadronic 
effects. We can answer to this question by performing a 
SM fit of the relevant experimental data (i.e. the values of 
ek, A mB d , V u b/V c b and the bound on Am Bs /Am Bd ) with 
ek itself excluded from the data to be fitted. It has been 
recently noticed [49, 50] that this kind of fit gives an inter¬ 
esting result: even omitting ek (the only so far observed 
CP-violating effect) a ‘good’ fit is possible only in presence 
of CP violation. 

We answer to the question in fig. 7, where we show 
the best-fit values of the Wolfenstein parameters (p,g). In 
the improved Wolfenstein parameterization of the CKM 
matrix, 


Kb = | Kb | e" i7CKM =A\\p- ig) 

V td = \V td \e~ i0CKM =A\ 3 (l-I\ 2 )(l-p-ig). 

In table 2 we list the values of the parameters used in the 
fit (we use the same notations and values of [50]). If we 
treat the errors on these parameters as standard deviations 
of Gaussian distributions, we obtain the fit shown in fig. 7a 
(ek included in the fit) and 7b ( ek not included in the fit). 
If we instead assume that the parameters p with theoretical 
uncertainty A p have a flat distribution with same variance 
as the gaussian one, the result of the fit is essentially the 
same. The case g = 0 (no CP violation in the CKM matrix) 
fits the data worse than other values, but is not completely 
outside the 99% ‘best fit’ region, as found in [50]. 

Since there are only few data to be fitted and the ‘good 
fit’ region is not very small it is not completely safe to use 
the standard approximate analytic fitting tecniques [51] 
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Figure 8: Naturalness distribution of various possibly interesting ‘non minimal’ supersymmetric effects (in the 
upper row, from left to right: BR(/j —► eff), \d e \/(e • cm), | = | arg(Ams d /Am^ I )|, direct CP asymmetry 

in b —> sj; in the lower row: supersymmetric correction to the CP asymmetry in the Bd —» 4>Ks decay, to 
dw /(e • cm), to ek, to the mixing induced CP asymmetry in b —► s'y) in unified supergravity with rjf = 0.5 (all 
assumptions are listed in the text). The vertical axes contain the values of the loop effects. The light gray part 
of the distributions comes from spectra where one of the loop effects is too large, the medium gray part from 
spectra accompanied by a discoverable effect; the black part from spectra where all the loop effects are too small. 
The horizontal dashed lines in each plot delimit the smallest effect that we estimate detectable (the continuous 
lines delimit already excluded affects). 


(we have performed our fit using a Monte Carlo technique). 
More importantly, in such a situation the exact result of 
the fit in general depends on the choice of the parameters 
to be fitted — for example the CKM angles /3ckm , 7ckm 
instead of p , g. This dependence becomes more important 
when studying if the experimental data on A ms d , A raj, 
and |Vuf,/F c (,| allow the CKM matrix to be real (77 = 0). 
To overcome all these problems we directly study how well 
any given particular value of (p, 77 ) is compatible with the 
experimental data, irrespective that other values could fit 
the data better or worse. When there are few experimental 
data this question is not necessarily numerically equivalent 
to asking what values of the parameters give the ‘best fit’. 

Again the results of this kind of analysis do not de¬ 
pend significantly on the shape (Gaussian or flat) of the 
distribution of the parameters with dominant theoretical 
error. The result is shown in fig. 7c, assuming Gaussian 
distributions for all uncertainties. We see that if 77 = 0 and 
p ~ 0.3 there is no unacceptable discrepancy between the 
experimental data and the theoretical predictions (using 
flat distributions 77 = 0 and p ~ 0.3 would be perfectly 
allowed). We conclude* that we cannot exclude a SM con¬ 
tribution to £k in the range (—2.5-v2.5)e“ p . Consequently 
we will consider allowed a supersymmetric correction to ek 
smaller than 3.5e^ p . 


*We however mention that a preliminary Tevatron study of 
the CP asymmetry in the decay Bd —> ipKs [52] disfavours 
negative values of 77 . 


Detectable supersymmetric effects 

The next question is: how sensitive to new physics will be 
the experiments performed in the near future? A planned 
experiment at PSI is expected to explore BR(p —> eff) 
down to 10 ~ 14 [53]. It seems possible to improve the search 
for the electron EDM (and maybe also the search for the 
neutron EDM) by an order of magnitude [54]. 

Concerning B and K physics, various new experiments 
will be able to measure CP asymmetries accurately. With 
the possible exception of A^f 7 that is difficult to mea¬ 
sure, the discovery potential is however very limited by the 
fact that the SM background has large QCD uncertainties. 
For example the precise measurement of the phase of B 
mixing (obtainable from the CP asymmetry in Bd{t) —> 
ifKs) PB d = 2/3ckm + <Pfld SY , says neither the value of 
/3 ckm = argVfd nor if a supersymmetric effect is present, 
0 - All proposed strategies that allow to disentan¬ 
gle the supersymmetric contribution from the SM back¬ 
ground suffer of disappointingly large (few 10 %) hadronic 
uncertainties: 

• One way for searching a supersymmetric effect in B 
mixing is the following. If one assumes an approxi¬ 
mate SU (3) symmetry between the three light quarks 
(u, d, s ) it is possible to reconstruct the unknown SM 
gluonic penguins (that affect the decay modes that 
allow to separate /3 ckm and <p^f d SY ) from the branch¬ 
ing ratios of the decays Bf —* n + K°, B d —> n~K + , 
B° d -► ^+7T- and their CP conjugates [55]. In this 
way it seems possible to detect a 30% correction to 
the phase of B^B^f mixing [39]. 

• Another way for searching a supersymmetric effect 
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in B mixing is the following. In the SM the di-lepton 
asymmetry in BdB d decays [56, 57] is given by 


T-'biV 

A^ = Im Fl2 


Mf 2 M 


0.001 arg 


Mf 2 M 


10 “ 


(M 12 — iF 12/2 is the off-diagonal element of the ef¬ 
fective Hamiltonian in the ( B , B ) basis) and is sup¬ 
pressed due to a cancellation between contributions 
of the u and c quarks. This cancellation could be sig¬ 
nificantly upset by unknown QCD corrections [56]. 
If this does not happen lepton asymmetries are a 
useful probe for a SUSY effect [57]: 


pSM -| n -3 

AjV SY = Im „ a „ ri2 „ gtIOV « 


• From a SM fit of the future precise experimental data 
(i.e. using, for example, a measurement of Amj, and 
of the CP asymmetry in Bd —> %j)Ks) it will be possi¬ 
ble to predict the SM value of ek with a uncertainty 
of about 25%, mainly due the future uncertainties 
on the Wolfenstein parameter A and on the matrix 
element of the operator that gives ek in the SM, 
determined with lattice techniques. A real improve¬ 
ment in the lattice computation will come only when 
it will possible to avoid the ‘quenching’ approxima¬ 
tion. We estimate (maybe a bit optimistically) that 
it will possible to detect supersymmetric corrections 
to ek larger than 30% • £™ p . It will also be possible 
to try to detect a supersymmetric corrections to the 
phase of B mixing from a global fit of the future ex¬ 
perimental data; it is not possible now to say if this 
technique can be more efficient than the direct ones 
discussed above. 

A theoretically clean way of detecting a supersymmetric 
correction to CP violation in BB mixing would result from 
a precise determination of the phase of B s mixing (that is 
very small in the SM; but its measurement does not seem 
experimentally feasible) or of BR(/7 + —> n + vv)/BR{K l —> 


parameter 

value 

Am Bd 

Am b s 

£K 

v ub /v cb 

m t (m t ) 

(0.471 ±0.016)/ps 
> 12.4/ps, see [37] 
(2.28 ±0.02) • 10~ 3 
0.093 ±0.016 
(166.8 ± 5.3) GeV 

m c (m c ) 

A 

Bk 

BTfB 

i 

hi 

V2 

773 

r/B 

(1.25 ±0.15) GeV 
0.819 ±0.035 

0.87 ±0.14 

0.201 ± 0.042 

1.14 ±0.08 

1.38 ±0.53 

0.574 ± 0.004 

0.47 ±0.04 

0.55 ±0.01 


Table 2: Values of parameters used in the determina¬ 
tion of (p, rf). The parameters in the upper rows have 
mainly experimental errors; the parameters below the 
middle horizontal line have dominant theorethical er¬ 
rors. 


n°uv) (that has small QCD uncertainties [58]). Hopefully 
these decay rates will be measured with ~ ± 10 % accuracy 
in the year 2005 [59]. Apart for this possibility we do not 
know any way that allows to detect corrections to B and 
K mixing from new physics smaller than (20 -j- 30)%. As 
a consequence we cannot be sure that a precise measure¬ 
ment of the CP asymmetry in B —> 1 pKs really measures 
the CKM angle /3ckm or if instead it is contamined by a 
~ 10 % new physics contribution. 

To summarize this long discussion, we have added to 
each plot of fig. 8 a horizontal dashed lines that delimits 
the smallest effect that we estimate detectable (while the 
continuous lines delimit already excluded effects). 

Discussion 

We have computed the loop effects characteristic of ‘non 
minimal’ supersymmetry assuming a 12/3 splitting between 
the three generations of sfermions with 77 / = 1/2 (see 
eq. (5)). We now discuss what happens if we modify our 
assumptions. 

The minimal effect that motivates 77 / 1 could give a 

much larger effect, 77 -C 1. This limiting case is interesting 
also for different reasons [36]. If 77 / « 0 leptonic effects 
are too large in about 80% of the parameter space and 
are almost always discoverable by planned experiments. 
Among the hadronic observables, a detectable supersym¬ 
metric effect is sometimes contained in ek and djv, while 
CP violation in B mixing and decays ( b —> sy, b —> sss) is 
interesting only for some values of the parameter that pro¬ 
duce too large effects in the other leptonic and hadronic 
observables. 

Alternatively, the 12/3 mass nondegeneration could be 
smaller. If 77 = 0.9 the GIM-like cancellation is so strong 
that no effect exceeds the experimental bounds, but lep¬ 
tonic effects remain discoverable in almost 50% of the pa¬ 
rameter space. Some hadronic effects still have a small 
possibility of being discoverable ( ek , dw, CP violation in 
B decays). If 77 = 0.95 only the leptonic signals have a low 
probability (~ 5%) of being discoverable. 

Apart for the value of 77 , our simplifying but question¬ 
able assumptions can be incorrect in different ways: 

1. We have assumed that all sfermions of third gen¬ 
eration are lighter than the corresponding ones of 
the first two generations. Maybe only some type of 
sfermions have a significant 12/3 non-degeneration: 
for example the ones unified in a 10 -dimensional rep¬ 
resentation of SU(5). In this case BR(/x — * ey) gets 
suppressed by a factor ~ (m p /m,) 2 , but likely re¬ 
mains the only interesting signal of new physics. It 
could instead happen that the squarks (but not the 
sleptons) have a significant 12/3 non-degeneration. 
In this unmotivated case a detectable effect in B 
physics is possible and accompanied always by a de¬ 
tectable correction to ek and often by a neutron 
EDM larger than 10 -26 e ■ cm. 

2. We have assumed fermion/sfermion mixing angles of 
the same order as the CKM ones. Recent neutrino 
data [60] suggest the presence of a large mixing angle 
in the leptonic sector. It is easy to compute how 
the various effects get enhanced in presence of some 
mixing angle larger than the CKM ones assumed in 
eq.s ( 6 ). 

3. We have assumed that the A terms alone do not in¬ 
duce interesting effects. In unification models the 
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masses of light quarks and leptons do not obey sim¬ 
ple unification relations: this suggests that they arise 
from higher dimensional operators with complex gauge 
structure. In this case, the A terms cannot be univer¬ 
sal and can induce a new large effect in /i —> ey [61] 
and a maybe detectable CP asymmetry in b —> 57 . 

In none of these cases CP violation in B mixing is an in¬ 
teresting effect. Even if supersymmetry gives rise to larger 
effects than the motivated ones that we have here studied 
an effect in B mixing is severely limited by the necessity of 
avoiding too large effects in leptonic observables, in EDMs 
and in ek (but cannot be excluded, because all bounds can 
be avoided in one particular situation). 

5 Conclusions 

In conclusion, the negative results of the recent searches 
for supersymmetric particles, done mostly at LEP2, pose a 
serious naturalness problem to all ‘conventional’ supersym¬ 
metric models. Fig. 1 illustrates the problem in a simple 
case where it is as mild as possible. Why the numerical val¬ 
ues of the supersymmetric parameters should he very close 
to the limiting value where electroweak symmetry breaking 
is not possible? 

There are two opposite attitudes with respect to the 
problem, that give rise to different interesting conclusions. 

One may think that supersymmetry has escaped detec¬ 
tion due to an unlucky numerical accident. This happens 
with a ~ 5% probability (or less in various particular mod¬ 
els), so that this unlucky event is not very unprobable. If 
this is the case we can study the naturalness probability 
distribution of supersymmetric masses in the small remain¬ 
ing allowed range of parameters of each model. It is no 
longer very unlikely that the coloured sparticles have mass 
of few TeV due to an accidental cancellation stronger than 
the ‘minimal one’ necessary to explain experiments. A sec¬ 
ond accident — just as unprobable as the one that has pre¬ 
vented the discovery of supersymmetry at LEP2 — could 
make the same job at LHC (assuming that it will not be 
possible to detect coloured sparticles heavier than 2 TeV. 
When the coloured sparticles are so heavy, the charginos 
and neutralinos are also too heavy for being detectable via 
pair production followed by decay into three leptons [62]). 
Even so, LHC has very favourable odds of discovering su¬ 
persymmetry. Before LHC, we estimate that there is a 
40% probability of detecting a supersymmetric correction 
to BR(B —> X s -y) and a 10% probability that supersymme¬ 
try affects the anomalous magnetic moment of the muon 
in a detectable way. Other interesting supersymmetric sig¬ 
nals are naturally possible only if supersymmetry gives rise 
to new flavour and CP violating phenomena. A detectable 
effect of this kind can be present almost everywhere for 
appropriate values of the hundreds of unknown relevant 
supersymmetric parameters. For this reason we have con¬ 
centrated our attention to a subset of strongly motivated, 
and thus controllable, signals [15]. We find that fi —> ey 
and the electron EDMs are interesting candidates for a 
supersymmetric effect. Effects in the hadronic sector are 
possible but not very promising (the neutron EDM and ek 
seem more interesting than CP violation in B mixing and 
decays: see fig.s 8 ). 

On the other hand one may instead think that 5% is a 
small probability: after all 95% is often used as a confidence 
probability level for excluding unseen effects. If one still be¬ 
lieves that supersymmetry at weak scale solves the SM nat¬ 
uralness problem, the supersymmetric naturalness problem 


motivates the search of unconventional models that natu¬ 
rally account for the negative results of the experimental 
searches. The problem would be alleviated by an appropri¬ 
ate correlation, for example between fi and M 3 . We know 
no model that makes this prediction; the naturalness prob¬ 
lem implies that any model really able of predicting /x/A /3 
has a large probability of making a wrong prediction. Even 
if some high energy model predicts the desired cancellation, 
this delicate cancellation will not survive at low energy due 
to large RGE corrections. In other words the required value 
of /X/M 3 depends on the values of 0 : 3 , At,... (for example 
we need /x(Mgut)/M 3 (Mgut) ~ 1.5 if AtG = 0.5, and 
« 2.5 if Ate = 1). 

Models where supersymmetry is mediated at lower en¬ 
ergy can thus have some chance of being more natural. 
However, exactly the opposite happens in the only appeal¬ 
ing models of this kind: gauge mediation models. These 
models predict that the right handed sleptons are lighter 
than the mass scale in the liiggs potential, so that natural¬ 
ness problem is stronger than in supergravity models. More 
importantly, if the supersymmetry breaking scale is so low 
that the neutralino decays in a detectable way, both LEP 
and Tevatron experiments give so stringent experimental 
bounds on the neutralino mass that make this scenario un¬ 
natural. 

A different way of alleviating the problem consists in 
having a sparticle spectrum more degenerate than the ‘con¬ 
ventional’ one. Since the Tevatron direct bound on the 
gluino mass is weaker than the indirect bound obtained 
from LEP2 assuming gaugino mass universality, it is pos¬ 
sible to reduce the mass of coloured particles (and conse¬ 
quently their large RGE corrections to the Z mass, that 
sharpen the naturalness problem) by assuming that gaug¬ 
ino mass universality is strongly broken. This possibility 
have recently been discussed in [63] in the context of su¬ 
pergravity. A gauge mediation model that gives a sparti¬ 
cle spectrum different from the conventional one (by as¬ 
suming an unconventional messenger spectrum) has been 
constructed in [64], In both cases, the more degenerate 
sparticle spectrum allows to reduce significantly the orig¬ 
inal FT parameter. Still, we believe that these solutions 
do not completely remove the unnaturalness. Even if all 
the sparticles and the W and Z bosons have now arbitrary 
but comparable masses, why only the SM vector bosons 
and not one of the many (~ 10 ) detectable sparticles have 
been observed with mass below Mz? It will be possible to 
concretely explore this possibility at Tevatron in the next 
years. 

Finally, a more original (but apparently problematic) 
approach is discussed in [65]. To conclude, we know no 
model that really predicts that sparticles are heavier than 
the Z boson (while we know many models that make the 
opposite prediction). 
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A Experimental bounds 

We now summarize the present experimental bounds on su¬ 
persymmetric particle masses. The informations are mainly 
extracted from [ 1 , 2 ], 


where 

V = (a 3 (nw) /«3(Ms)) 1/2b5 

and 65 = 23/3 is the coefficient of the QCD (3 function at 
one loop with 5 flavours. The hadronic matrix elements of 
the operators are 


LEP2 bounds The most important bound is the one on 
the chargino mass: m x > 91 GeV unless M x — Mn < 4 GeV 
or mo < Mw- The bound on the higgs mass is also very 
important if tan /3 </ 3 is low and can be approximated with 

f 91 GeV if tan /3 < 4 
mh > y 83 GeV if tan f3 > 4 

The bound on charged sleptons is > 80 GeV. 

Finally there is a very stringent bound on the neu- 
tralino mass: if the neutralino decays into 7 and grav- 
itino within the detector (like in gauge mediation models 
with Mm & 10' GeV and the neutralino is the ‘NLSP’) then 
mjv > 91 GeV. Tevatron experiments give a similar bound, 
Tn n > 75 GeV. 


Tevatron bounds The strongest bound on coloured 
sparticles come from Tevatron experiments. The bound on 
the gluino and squark masses can be approximated with 

/ 180 GeV if mq » M 3 
3 > \ 260 GeV if m<j « M 3 

The bound on the stop mass is much weaker: mi <; 75 GeV. 
In our analysis we have imposed that the pole mass of the 
top quark be in the range (175 ± 10) GeV. 


(- B d \0 LL \B£) 
{B d \ 0 RR \B 3 ) 


(B d \Ol R \B d ) 


{B d \Ol R \B d ) 


-^f\mBB L L 

-^S\"m, B B RR 




4 (m b ± m d ) 


\flm B Bl R \ 


4 (m b + m d ) 


2 


2 


for B d B d mixing. The corresponding expressions for B s Bs 
mixing can be obtained replacing d —► s. The correspond¬ 
ing expressions for KK mixing can be obtained replacing 
B —> K and b —> s. The parametrization of the matrix el¬ 
ements that we employ, different from the vacuum satura¬ 
tion approximation (VSA), is more convenient for a lattice 
computation of the B factors. Notice that the very large 
scale dependence arising from the factor r /~ 16 is largely can¬ 
celled by the scale dependence of the quark masses present 
in the hadronic matrix element. The enhancement of ef¬ 
fects in ek due to QCD evolution is thus less dramatic 
than argued in the first of ref. [43], as confirmed by the 
fact that all the i?_Lij(/xs)-factors turn out to be not much 
larger than 1 at fi B = mb{rrib)- 

The B factors computed in quenched approximation on 
the lattice are 


B CP violation in BB and KK 
mixing and supersymmetry 


We can write the effective Hamiltonian for A F = 2 (F = 
S, B) processes as 

H?/= 2 = CllO ll +C rr O rr +CI r OI r +CZ r OI r + h.c. 


where the relevant A F = 2 operators are 

Oll = (6L7M d i)(^n7M d i) 
O rr = (&k7M<4)(47^4) 

®lr = (&z,dk)(44) 
o x LR = (6l4)(44) 


Here i and j are colour indexes and for simplicity we have 
listed the operators relevant for B d B d mixing — the opera¬ 
tors relevant for B s Bs and KK can be obtained with triv¬ 
ial replacements of the quark flavours. The most general 
A F = 2 Hamiltonian contains two more operators that we 
have not considered because they are generated with neg¬ 
ligibly small coefficients in the supersymmetric scenario in 
which we are interested. 

The leading order QCD evolution of the coefficients 
C(/x) from nw = O(Mw) to Hb = 0(ms) in the SM 
is [43] 


Cll{17b) 

= ti^Cll^w) 

C RR ([1b) 

= V 4 C RR {l l w) 

Cl R (^b) 

= ri~ 16 Clr(hw 

^L R (^b) 

= V C lr {vw) 


B ll {iib) = 0.9 ±0.1, Bz r (ii b ) = ^(1.15 ±0.10), 

B rr {hb) = 0.9 ± 0 . 1 , B£r(a»b)»|( 1.15 ± 0 . 10 ), 

in the B d and B s systems [44], and 

B ll (ii b ) = 0.57 ±0.06, BZr(jj.b) = 1.01 ± 0.06 

B rr {ll b ) = 0.57 ±0.06, B rr (jj,b) = 0.78 ± 0.11 

in the X-system [45]. In the case of the B systems the 
VSA value of some Bl R coefficients is the number different 
from one that we have sorted out. All the B coefficients 
— including the ones for the K system — have been given 
renormalized at the scale fj, B = 4.2 GeV. At the same scale 
the relevant quark masses are [46] 

m b (fi B ) = (4.2 ±0.2) GeV 

m s {fj,B) = (0.105 ±0.017) GeV 

and 7 « 0.954 for /rw = mt(mi). Notice that the fact that 
Bl R ~ 1 in the K system means that the large enhance¬ 
ment {rriK/m s ) 2 of the matrix elements of the LR super- 
symmetric operators with respect to the matrix elements 
of SM operator Oll, predicted by the vacuum insertion 
approximation, is confirmed by lattice calculations. In our 
formulas we have omitted all details necessary for a NLO 
computation: NLO supersymmetric effects are only par¬ 
tially known and of the order of the uncertainty on the B 
parameters. 

With our normalization of the hadronic states B and K 
f K = 0.160 GeV, f B = (0.175 ± 0.35) GeV 
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and the measurable quantities A mg, Am* and Ek are 
given by 


Am s = 2{B\'H / t£ =2 \B) 

A m K = 2Re{K\Hts S=2 \K) 

i7r/4 lm(K\H^ s = 2 \K) 

£k = e - — - 

v2A?7iic 

The values of the coefficients in the supersymmetric sce¬ 
nario considered in section 4.2 are, for B&B j mixing: 


Cll | SM 
C^| SUSY 

Crr\ SVSY 


Cl 


I SUSY 


/~i X I SUSY 
r, B 


T/ 2 t/*2 a 

VtdVtb 
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• 2.52 
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«3 


d L b L b L b L g M 2 

y 2 . y* 2 ~ 03 x 

d Rb R b R b R g M 2 

}({r 3 - ri 2 }, {r 3 - r 12 }) 


= V 4 ~ h V^rV*r.V*r 


_ 2 
«3 


d _L b L V d R b R * b L b L Vb R b R 3 j\^2 

x{-4B^ + }({£ 3 - 42 }, {r 3 - n 2 }) 


= V,zV*i. V* 


«3 


d I#&L V d R b R Vb L b L V hR b R 9 

x{20B^ + }({4 — 4 . 2 }, {y 3 — ri 2 }) 


where, in order to avoid long expressions, we have intro¬ 
duced the following compact notations: 


4 = —— 

3 “ Ml ' 


tii = 


"‘12J, 


M| 


m, 

OR 

r 3 = 


M|’ 


ri2 = 


"‘ 12 * 


Ml 


and 

/({ai±a 2 }) = /(ai)±/(a 2 ), 


{/i±/ 2 }(a) = /i(a)±/ 2 (a) 
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